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We demonstrate that for two canonically conjugate operators q,p ,the global correlation (qp + 
pq) — 2(g) (p) , and the local correlations (q) (p) — (q) and (p) (q) — (p) can be measured exactly by Von 
Neumann-Arthurs-Kelly joint quadrature measurements . These correlations provide a sensitive 
experimental test of quantum phase space probabilities quite distinct from the probability densities 
of q,p. E.g. for EPR states, and entangled generalized coherent states, phase space probabilities 
which reproduce the correct position and momentum probability densities have to be modified to 
reproduce these correlations as well. 
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Von Neumann- Arthurs-Kelly Joint Measure- 
ments of Conjugate Variables in Quantum Me- 
chanics . Correlations between conjugate observables, 
being rather different from Bell type correlations [l[ 
among commuting observables, is a largely unexplored 
area with possible fundamental importance. We present 
a method for exact measurement of local and global cor- 
relations between conjugate observables in quantum me- 
chanics. We use the Von Neumann-Arthurs-Kelly Joint 
Measurements, realizable as heterodyne measurements 
in quantum optics. As a first application these corre- 
lations are used to experimentally test proposed phase 
space probabilities and to constrain costruction of such 
probabilities so that they reproduce not only the quan- 
tum probability densities of conjugate observables but 
also their correlation. 

Von Neumann Q not only proposed a model Hamil- 
tonian for accurate measurement of a single observable q 
but also noted that joint approximate measurements of 
canonically conjugate observables q,p were possible with 
accuracy limited by the uncertainty principle . Arthurs 
and Kelly generalized the Von Neumann Hamilto- 
nian to realize such joint measurements and deduce a 
measurement uncertainty relation which was soon proved 
with full generality by Arthurs and Goodman Q . Braun- 
stein, Caves and Milburn [f| analyzed the optimum ini- 
tial conditions and deduced that the optimum Arthurs- 
Kelly joint distribution for q, p is just the Husimi Q func- 
tion Stenholm Q proposed a realizable quantum 
optics Hamiltonian for such measurements. The vari- 
ous facets of the uncertainty relations so revealed have 
been reviewed by Busch, Heinonen and Lahti Q. The 
most satisfying thing is that Heterodyne measurements 
in quantum optics [9( can practically realize joint mea- 
surements of any pair of rotated conjugate quadratures, 
(aexp {—i0)/\/2 + h.c, — iaexp (—i0)/\/2 + h.c.) where 
q,p are given in terms of the photon annihilation and 
creation operators a, a', 

q={a + at)/V2,p = (a - a f )/(n/2). 



Arthurs-Kelly Results. Their idea is that the sys- 
tem (with position and momentum operators q,p) inter- 
acts with an apparatus which has two commuting ob- 
servables Xi,X2 and approximate values of system posi- 
tion and momentum are extracted from accurate joint 
observation of x\,X2- The Von Neumann- Arthurs-Kelly 
interaction during the time interval (to, to + T) is, 



H = K(qp! +PP2), 



(1) 



where p\,p2 are canonical conjugates of Xi,X2 respec- 
tively, the coupling K is large, and T is small, with 
KT = 1. During interaction time, H is so strong that 
the free Hamiltonians of the system and apparatus are 
neglected. Arthurs and Kelly start with the system- 
apparatus initial state, 



tp(q,Xi,x 2 ,to) = <P(q)xi(xi)X2(x2) 



(2) 



where, (j>(q) is the system state and the apparatus state 
is given by, 

Xi(x!) = n-^b- 1 / 2 cxp(-x 2 /(2b 2 )), (3) 

X2(x 2 ) = tt- 1 / 4 (26) 1 / 2 cxp(-2& 2 :rj), (4) 

and b/y/2 is the uncertainty of x\ in the initial apparatus 
state. They solve the Schrddingcr equation exactly and 
obtain the final joint probability density of the apparatus 
variables to be just the Husimi function [(|, 



where 



P(x U X 2 ) = \{4>b,x u x 2 \4>)\ K^), 



(5) 



(9) = {2irtf)- 1 '*exp(iqx 2 -(x 1 -qy/(4b')) (6) 

is a minimum uncertainty system state centred at q = 
x\,p = X2- Note that for any value of b, (xi) = (q), and 
(^2) = (p)j but the dispersions in X\,X2 are larger than 
those for the corresponding system variables q,p, 



{Ax 1 ) 2 = (Aq) 2 + ti > , (Ax 2 ) 2 = (A P ) S 
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they obey the "measurement or noise " uncertainty rela- 
tion , (units H = 1 ), 



AxiAz 2 > 1 



(8) 



obtained by varying b. Here the minimum uncertainty is 
twice the usual " preparation uncertainty" . Arthurs and 
Goodman Q gave a beautiful proof of this fundamental 
uncertainty relation, independent of any particular choice 
of the Hamiltonian. Further, for the x\ distribution to 
approximate q distribution closely, we need b <C Ag; for 
the x 2 distribution to approximate p distribution closely, 
we need b » (Ap)^ 1 . 

Pifa) = f P(x u x 2 )dx 2 = (2Tr)- 1 ' 2 b- 1 



dq\^(q)\ 2 c X p(~(x 1 ~q) 2 /(2b 2 )) ->6_x> |0(*i)| a , (9) 



P 2 {x 2 ) = J P{x 1 ,x 2 )dx 1 = {2Ti) 1 ' 2 b 
J dp\4>(p)\ 2 cx V (-(x 2 -p) 2 (2b 2 )) |0(x 2 )| 2 ,(lO) 

where 4>(p) denotes the Fourier transform of (f>(q). 

Exact Measurement of Quantum Correlations 
Between Conjugate Variables. The above equations 
show that the exact position and momentum probability 
densities of the system are recovered by the Arthurs-Kelly 
measurement in the limits b — ¥ and b — > oo respec- 
tively,!. e. in two experiments with very different initial 
apparatus states. It is a pleasant surprise that the joint 
measurement can nevertheless give local and global cor- 
relations between q,and p exactly. We define, 



(p)(q) 



(A(q)p + pA(q)) 
2(A(q)) ' 



<?)(P) 



(A(p)q + qA(p)) 
2<A(p)> 



(11) 

where (A) denotes the quantum expectation value of 
a self-adjoint operator A, and the projection operators 
A(q),A(p) are defined by, 



A(q) = \q)(q\, A(p) = \p)(p\. 



(12) 



For a pure state \cj>) 

(£>(?) = 

<9>(p) = 



we have the explicit expressions, 
Re(<f>*(q)(-i)d<l>(q)/dq) 



l#z)l 5 



Re($*(p)(i)dfo)/dp) 



(13) 



(14) 



We shall see that the local correlations (p)(q) — (p),and 
(q)(p) — (q) can be measured exactly for arbitrary q and 
p respectively for appropriate values of b. The global 
correlation (qp+pq) — 2(q) (p) is in fact exactly measurable 
for any value of b. 



For the Arthurs-Kelly measurement we define as for a 
classical distribution, 



(x 2 )a-k(xi) = / x 2 P(xi,x 2 )dx 2 /Pi(xi), 



(xi)a-k(x 2 ) = J xiP(xx,x 2 )dxi/P 2 (x 2 ), 



{x 1 x 2 ) A _ K 



X\X 2 P{x\ , x 2 )dx\dx 2 



(15) 



(16) 



(17) 



Substituting the value of P(x\,x 2 ), and doing the inte- 
gral over x 2 we obtain, 

J x 2 P(x l7 x 2 )dx 2 = (bV^)- 1 J dqdq'(t>{q)<t>*(q') 

(z 1 - g ) 2 + (x 1 -g') 2 v d%- g Q 
ex P ( 713. > 



= Re 



Ab 2 

dq , (xi ~q) 2 s v .sd(j){q) 

ex P( ^ ( 18 ) 



6V2tt 



dq 



where S(q — q') is the Dirac delta function. Similarly, we 
obtain, 



J xiP{x\ 1 x 2 )dxx = b^/2/n 



Re / dpexp(-2b z (x 2 - P y)(j)*(p)i 



dp 



(19) 



Taking the limits of b going to and oo yield respectively, 
(x 2 )a-k{xi) ^ 6 ^o (p){q = xi), (20) 



{xi) A -k{x 2 ) 



(q)(p = x 2 ). 



(21) 



Thus we have proved that the quantum position proba- 
bility density and the local correlation (p)(q) — (p) can 
be measured exactly with the initial condition b — > 0; 
the quantum momentum probability density and the lo- 
cal correlation (q) (p) — (q) can be measured exactly with 
the very different initial condition b — > oo. A similar 
calculation shows that for any value of 6, 



(2x 1 x 2 ) A -k = (qp+pq), 



(22) 



the global correlation is exactly measured in the Arthurs- 
Kelly (A-K) measurement. Thus, the A-K measurements 
with b — > and b — > oo equip us with exact probability 
densities of position and momentum as well as their exact 
local and global correlations. 

Experimental test of phase space probabilities 
by correlation measurements. We demonstrate that 
exact measurement of the correlations is a valuable tool 
to discriminate between various phase space probabil- 
ity densities which may give exactly the same posi- 
tion and momentum probability densities. The tremen- 
dous progress initiated by research on Bell inequalities 
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and quantum contextuality [l|, and their extension to 
phase space [l2| teaches us that in 2N dimensional phase 
space, a positive density can have a maximum of N + 1 
marginals reproducing quantum probability densities for 
arbitrary states. (E.g. for N = 2, probability densities 
of (gi, 52), (pi, 92), (pi,P2) can be reproduced.) Of course 
we know that all marginals of Wigncr's quasi-probability 
distribution (l3j agree with the corresponding quantum 
probablities for the state |(£).But we shall only cosidcr 
positive densities as candidates for a probability inter- 
pretation. Dc Broglie and Bohm proposed a posi- 
tive phase space density which reproduces the quantum 
position probability density but fails to agree with the 
quantum momentum probability density [llj. The most 
general positive densities with two marginals reproduc- 
ing quantum position and momentum probabilities 14| . 
and with N + 1 marginals reproducin g th e corresponding 
quantum probabilites are also known [12j . Roy and Singh 
[la ] built a new causal quantum mechanics symmetric in 
q,p in which the phase space density obeys positivity 
and the marginal conditions on momentum and position 
probabilities . For example for N = 1, the two densities 

p e (^)H#z)| 2 |<^)| 2 

S ( f P dp'\4>(p')\ 2 - r dq'\4>(eq' ,t)A , (23) 

where e = ±1 clearly reproduce the quantum position 
and momentum probabilities as marginals. 



Pe (q,p)dp = |#?)| 2 ; / Pe (q,p)dq = \4>(p)\ 2 . (24) 



To demonstrate the discriminatory power of the quantum 
correlation measurements we shall use them in several 
concrete examples to test these two phase space densities 
( for e = ±1 ), as well as the correlationless phase space 
density |0(g)| 2 |0(p)| 2 , all of which reproduce quantum 
q,p probability densities. 

(i) Free particle spreading wave packets for non- 
rclativistic particle of mass m. At the time to of the 
A-K measurement, let 



<fi(p) = (ira) 1 ^ 4 exp[ 



(P-P) 2 
= f,(Ag) 2 = 



P 



2a "2m' 
1 + (at /m) 2 
2a ' 



(25) 



The Roy-Singh q,p symmetric causal quantum mechanics 
gives, for e = ±1, 

A» 

(p)(q)± - (P) = ±-rH«Z - Ph/m), (26) 
Ag 

whereas the Arthurs-Kelly correlation is, 



(X2)(Xl)A-K - (P) = , A ro— -To (Xl - Pto/m) 



(Aq) 2 + b 2 



and it's limit — >■ is the true quantum correlation. The 
correlationless phase space density gives for the above 
correlation and the e = — 1 case gives a negative corre- 
lation, both disagreeing with the quantum correlation, 
whereas the ratio of the correlation in the A-K measure- 
ment to that in the e = 1 Roy-Singh causal density (see 
figure) approaches unity for b/ Aq <C 1 and AqAp ^> 1/2; 
i.e. there is agreement with quantum mechanics only 
when the uncertainty product is large. Similarly, for cor- 




0.4 Aq 



FIG. 1: For the free particle expanding Gaussian wave packet, 
the ratio of the correlation (x2)(xi) — (x2) in the Arthurs-Kelly 
measurement to (p)(q) — (p) in the e = 1 Roy-Singh causal 
phase space density is plotted for various values of b/Aq and 
of AqAp.The causal correlation agrees with the quantum cor- 
relation (i.e. the b — > limit of the A-K correlation )only for 
large values of the uncertainty product. A convex combina- 
tion of the e = 1 and e = — 1 causal phase space densities can 
reproduce quantum correlations exactly. (Figure computed by 
Arunabha S. Roy.) 



relations at given p, the Roy-Singh causal quantum me- 
chanics gives, for e = ±1, 

( q )(p) ± -(q)=±^L(p-f3), (28) 

which agrees only for e = 1 and only for large AqAp with 
the quantum correlation which is the b — > 00 limit of the 
Arthurs-Kelly correlation , 



(X1)(X 2 )A-K - (q) = ( Ap )2 + ( 4 fr2)-l ( X2 _ Pi' ( 29 ) 

For the global correlation, the Roy-Singh causal quantum 
mechanics with e = ± gives, 



(2qp) ± -2(q)(p) =±2 AqAp, 



(30) 



of which only the e = 1 correlation agrees with the quan- 
tum correlation, 



(27) 



(qp + pq) - 2(g) (p) = ^(2AqApf - 1, (31) 
provided that 2 AqAp ^> 1 . 
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(ii) Generalized coherent states of light [Ta]. These are 
displaced n-th excited states of the oscillator of frequency 

w, 

1 _ p 

<t>n,a{(l,to) = (q-q\n) exp (-iut (n + -) +ip(q - -), 

(32) 

where a = j4exp (— i(wio + 9)), q = Rea,p = Ima, and 
A, 9 are real constants. Here the Roy-Singh causal quan- 
tum mechanics with e = ±1 gives, 

(p)(q)±-(p)=±(q-q), (q)(p)± - (q) = ±(p - p), (33) 



(2qp) ± -2(q)(p)=±(2n + l). 



(34) 



In contrast quantum mechanics gives zero for the above 
three correlations and thus agrees with the correlationless 
phase space density. 

Phase space probabilities reproducing quantum 
position and momentum probabilities and corre- 
lations exactly. Surprisingly, in both the examples 
considered above, convex combinations of the Roy-Singh 
phase space densities with e = ±1, 

p(q,p)c = Kp+{q,p) + *-p-{q,p), (35) 



< A± < 1, A+ + A_ = 1. (36) 

where the state dependent constants A± are chosen to re- 
produce the quantum global correlation (qp+pq)~2(q) (p) 
yield local correlations also equal to the corresponding 
quantum local correlations. Explicitly, in cases (i) and 
(ii) of Gaussian packets and generalized coherent states, 

1 1 



(i) A± = - ± - y/i - (2A 9 Ap)- 2 , (ii) A± = 1/2. (37) 

EPR states. A normalizable version of the original 
EPR state [13] \qi~ 92 = qo)\pi+P2 = Po) of two particles 



<i>{qi - qi-.pi +P2) = <Pi(qi - 92)<febi +P2), (38) 

where the individual Gaussian wave functions, 

Mn-<h) = (^r^cxpe- ^-f ~ go)2 ), (39) 



1 i , ^ 1 ^-l/4 t (P1+P2-P0), , An . 

<h(Pi +P2) = (na. 2 ) 1 exp( ) (40) 

2a 2 

are sharply peaked at q\ — 92 — qo = and p\ +p 2 — Po = 
respectively in the limits a\ — > 0, «2 — > 0. We now 
construct the phase space density, 

, (pi -P2), Aqi + q2) , . „ 
pw(qi-q2, s )P2C{ o ,(pi+P2)), (41) 



with the two factors p\c and p 2 c made to fit respec- 
tively the (qi -q 2 , (pi ~P 2 )/2) and (91 + q 2 )/2, (p 1 +p 2 )) 
correlations in the Gaussian states 4>i , <p 2 using convex 
combinations of the Roy-Singh phase space densities de- 
scribed above. This phase space density reproduces ex- 
actly, the above quantum correlations as well as quantum 
joint probability densities of the four commuting pairs of 
variables 91-92, (qi+q 2 )/2; 91-92,^1+^2; {ox+Qa) /^,Pi~ 
P2\P\ +P2, (Pi -P2)/2. 

Entangled generalized coherent states. For two 

modes of light with the same frequency an exact solution 
of the Schrodinger equation at time to is the entangled 
generalized coherent state , 



0m,a((?i + «a)/V2,to)^((?i-ffi!)/A*o). (42) 



where m, n are integers, a, j3 complex constants and the 
factors m . Q , <t>n,p are generalized coherent states defined 
before . A phase space probability reproducing the rel- 
evant quantum correlations and probabilities exactly is 



,(31 + 92) (P1+P2K ,(91-92) (pi 

PmCl 7= , 7= )PnC{ 



Pi) 



V2 



V2 



V2 



V2 



(43) 

where p m c , PnC are arithmetic means of the e = ± Roy- 
Singh phase space densities for the states <p m .a,<Pn.p re- 
spectively. 

Future directions. The central point is the exact 
measurability of local and global correlations between 
conjugate observables. Actual joint quadrature measure- 
ments to test their correlations will be very interesting. 
The Arthurs-Kelly joint measurements and hence the 
possibilities of exact measurements of quantum correla- 
tions between conjugate variables can be generalized to 
2iV-dimcnsional phase space. An interesting question is 
triggered by the success in exact reproduction of cho- 
sen quantum correlations and probabilities in the special 
states (including entangled states) considered. Can we 
construct phase space probabilities reproducing quantum 
position and momentum probabilities and their correla- 
tions exactly for every quantum state ? 

I wish to thank Arunabha S. Roy for the crucial fig- 
ure on experimental test of correlations , Virendra Singh 
for many intensive discussions and the Indian National 
Science Academy for financial support. 
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